Time scales in the dynamics of an interacting quantum dot 
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We analyze the dynamics of a single-level quantum dot with Coulomb interaction, weakly tunnel 
coupled to an electronic reservoir, after it has been brought out of equilibrium, e.g. by a step-pulse 
potential. We investigate the exponential decay towards the equilibrium state, which is governed 
by three time scales. In addition to the charge and spin relaxation time there is a third time scale 
which is independent of the level position and the Coulomb interaction. This time scale emerges in 
the time evolution of physical quantities sensitive to two-particle processes. 



PACS numbers: 73.23.-b,73.23.Hk,73.63.Kv 

I. INTRODUCTION 

The control and manipulation of single electrons in 
mesoscopic systems constitutes one of the key ingredi- 
ents in nanoelectronics. The study of single-electron 
sources^— in the high-frequency regime has attracted 
a great interest due to their potential application in 
quantum electron optics experiments, in metrology, and 
in quantum information processing based on fermionic 
systems^— In this work we study the time evolu- 
tion of a quantum dot (QD) tunnel coupled to a sin- 
gle electronic reservoir, as depicted schematically in 
Fig. [TJa). In the presence of some time-dependent 
voltage modulations, this system defines the building 
block of the typical single-electron source, namely the 
mesoscopic capacitor^ In the linear-response regime, 
the relaxation behavior of such a mesoscopic capaci- 
tor has been extensively studied theoretically^^— and 
experimentally^ revealing the quantization of the charge 
relaxation resistance! 16 i 17 ' 19 i 21 : 22 i 25 ' 26 On the other hand, 
the application of nonlinear periodic potentials to the 
mesoscopic capacitor yields the controlled emission and 
absorption of electrons at giga-hertz frequencies From 
these experiments the average charge as well as current 
correlation o 7 ' 27 ' 28 after each cycle of the potential ap- 
plied have been extracted. These results demonstrate 
the importance of investigating the dynamics of this 
kind of single-electron sources. In some of the recent 
realization a 2 ' 7 ' 26 the Coulomb interaction is weak; how- 
ever, in small-sized QDs the Coulomb blockade is, in gen- 
eral, strong and it is, therefore, desirable to include it in 
the theoretical analysis ^ 25 ' 29 since it may even domi- 
nate time-dependent phenomena, see e.g. Ref. [30L The 
time-evolution of interacting quantum dots after the cou- 
pling to the leads has been switched on, has, e.g., been 
studied in Refs. I31U34I and references therein. 

Here we investigate the exponential relaxation of a QD 
towards its equilibrium state after its has been brought 
out of equilibrium by applying, e.g., a voltage step pulse. 



We consider a voltage pulse that affects the occupation 
of only a single orbital energy level. The level can be spin 
split due to Coulomb interaction. In an earlier work^ 
some of the present authors investigated the decay of 
charge and spin of such a single level QD. It was found 
that the relaxation of charge and spin are given by rates 
which differ from each other due to Coulomb repulsion. 
Since the reduced density matrix of a QD with a single 
orbital level with spin is four dimensional, there are thus 
three rates which govern the relaxation of the diagonal 
elements of the density matrix towards equilibrium (plus 
one which is always zero and corresponds to the stable 
stationary state). In addition to the rates that govern 
charge and spin there is a third rate that appears in the 
relaxation of a single level QD with spin and with inter- 
action. This additional rate is the subject of this paper. 

Interestingly, this additional time scale is independent 
of the interaction and of the dot's level position. It is 
shown to be related to two-particle effects and appears, 




FIG. 1: (Color online) Schematics of the models: a) Single 
level QD with Coulomb interaction, U, coupled to a normal 
lead with a tunneling strength F. Dot occupations can be 
measured via the current passing through a nearby quantum 
point contact (QPC) capacitively coupled to the dot. b) QD 
attached to an additional superconducting contact, c) QD 
coupled to a ferromagnetic lead. 
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e.g., in the time-evolution of the mean squared devia- 
tions of the charge from its equilibrium value. We study 
in detail the perturbations leading to a relaxation of the 
system with the additional decay rate only, and find that 
it is indeed related to two-particle correlations. We also 
propose a procedure to separately read out the different 
relaxation rates occurring in the dynamics of the QD ex- 
ploiting the sensitivity of a nearby quantum point contact 
to the occupation of the QD, sec Fig. [1] (a). 

In order to further clarify the properties of the ad- 
ditional time scale, we extend our study to two other 
setups: a QD proximized by an extra, superconducting 
electrode and tunnel coupled to a normal lead; and a QD 
tunnel coupled to a ferromagnetic lead, see Fig. Q] (b) and 
(c). 



II. MODEL 

We consider a quantum dot coupled to an electronic 
reservoir. We assume that the single-particle level spac- 
ing in the dot is larger than all other energy scales, so 
that only one, spin-degenerate level of the QD spectrum 
is accessible. At a certain time to the system is brought 
out of equilibrium, e.g. by applying a gate potential, 
and afterwards relaxes to an equilibrium dictated by the 
Hamiltonian H = ifrj + Ht + H rcs . The Hamiltonian H-q 
of the decoupled dot 



(1) 



contains the spin-degenerate level e and the on-site 
Coulomb energy U for double occupation of the dot. The 
creation (annihilation) operator of an electron with spin 
(7 =f,l on the dot is denoted by d) a (d a ) and h a is the 
corresponding number operator. The reservoir is mod- 
eled by the Hamiltonian -ff rcs = e kc\ cr c ka , in which 

ct (c ka ) creates (annihilates) an electron with spin a and 
momentum k in the lead. The coupling between the dot 
and the reservoir is described by the tunneling Hamilto- 
nian Ht = J2 k a {Vc k(T d a +H.c), where V is a tunneling 
amplitude, which we assume to be independent of mo- 
mentum and spin. By considering a constant density of 
states v in the reservoir, the tunnel coupling strength T 
is defined as T = 2itv\V\ 2 . 

In the remainder of this paper, we focus on the re- 
laxation behavior of the quantum dot to its equilibrium 
state and in particular on how this relaxation manifests 
itself in measurable quantities. We are not interested 
in the dynamics of the reservoir, thus the trace over its 
degrees of freedom is performed to obtain the reduced 
density matrix of the QD. The Hilbert space is spanned 
by the four eigenstates of the decoupled dot Hamiltonian, 
{|x)}, where |0) represents the unoccupied dot, the dot 
is in the state \a) when being singly occupied with spin 
(7 =t,|, and \d) is the state of double occupation. The 
energies related to these states are Eq = 0, E a = e and 



Ed = 2e + U, where we set the electrochemical potential 
of the reservoir to zero. As we consider spin-conserving 
tunneling events, the off-diagonal elements of the reduced 
density matrix evolve independently of the diagonal ones 
(which are the occupation probabilities). We can, there- 
fore, consider these probabilities alone, which arranged 
in a vector are given by P = (po,Pf,Pi,Pd) T and fulfill 
the condition J2jPj(t) = 1- The time evolution of the 
occupation probabilities is governed by the generalized 
master equation 



dP(t) 
dt 



W{t,t')P(t')dt' , 



(2) 



In 



where the matrix elements W X:X '(t,t') of the kernel 
W(i, t') describe transitions from the state |x') at time 
t' to a state \\) at time t. 

We consider now the dynamics of the system after be- 
ing brought out of equilibrium at time to- Since for t > to 
the total Hamiltonian is time independent, the transition 
matrix elements depend only on the time difference t — t', 
i.e. W(t, t') — > W(t— t'). Furthermore, we are interested 
in the exponential decay towards equilibrium. To be 
more specific, we will therefore consider only the leading, 
time- independent, prcfactor of the exponential functions. 
Time-dependent corrections to the pre-exponcntial func- 
tions, that generally may appear ; 34 ' 35 are disregarded. 
Furthermore, when focussing on times t distant from the 
switching time to, such that the difference t — to is hence 
much larger than the decay time of the kernel W(t — t'), 
we can replace the lower limit of the integral in Eq. ([2]) by 
— oo. Expanding the probability vector P(t') in Eq. @ 
around the measuring time t we find 2 ^ 



dt ^-i n\ dt n 



(3) 



n=0 



Here we introduced the Laplace transform of the kernel 

w(z) = f ^ w(t - t'y-^-^dt 1 ', with w = w(z)| z=0 

and the n-th derivative of the kernel with respect to the 
Laplace variable 9™W = [d n 'W(z)/dz n ] z _ . The formal 
solution of Eq. ([3]) is given by 

P(i) - cxp(At)P in , (4) 

which depends on the initial probability vector P ln at 
t = to, where the initial values for the system parameters 
are given by the ones just after the switching time t . 
The matrix A includes Markovian and non-Markovian 
processes^ In the following, we consider the limit of 
weak coupling between quantum dot and reservoir and 
limit ourselves to a perturbation expansion up to sec- 
ond order in T, which is valid for the regime where the 
tunnel coupling T is much smaller than the energy scale 
set by the temperature ksT. The perturbativc expan- 
sion of A is A = AW + A< 2 ) with A« = WW and 
A( 2 > = W^+clW* 1 ) - WW, where the number in the su- 
perscript represents the power of T included in the transi- 
tion matrix W. Notice that the first non-Markovian cor- 
rection, i.e. the term dW^ ■ W^ 1 ' is present in second- 
order in the tunnel coupling. The evaluation of the kernel 
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within a perturbativc expansion can be performed using 
a real-time diag rammatic technique ) 37,38 which has been 
used in Ref. [29| in order to extract the exponential decay 
of spin and charge in the system studied here. Consid- 
ering Eq. Q , we see that the rates defining the decay of 
the state into equilibrium are found from the eigenvalues 
of the matrix A, which turn out to be real and non- 
positive. The matrix A is not Hermitian, as expected 
since we deal with a dissipative system, and hence has 
different left and right eigenvectors, 1^ and rv 



III. RESULTS 

A. Relaxation of the reduced density matrix 

The time-dependent probability vector, P(£), can be 
expressed in terms of the right eigenvectors of A, each 
being related to a decay with a different rate. The left 
eigenvectors determine the observable that decay with a 
single time scale only, see also the appendix. 
In the following we discuss the exponential relaxation 
towards equilibrium of the vector of occupation proba- 
bilities, in first order in the tunneling strength T. 



1. Nomnteracting limit 

We start by briefly discussing the simplest case of a 
single spinless particle. This limit is obtained, when 
a magnetic field much larger than the temperature is 
applied, B k^T. The Hilbcrt space of the system 
is two dimensional and spanned by the states |0) and 
|1) for the empty and singly-occupied dot respectively, 
whose occupation probabilities are arranged in the vec- 
tor P = (po,Pi) T - The decay to the stationary state 
is governed by matrix A^ (defined equivalently to A^ 1 ) 
but for the two-dimensional Hilbert space for the problem 
at hand) which contains a single relaxation rate, namely 
the tunnel coupling T, as intuitively expected. 

We now include the spin degree of freedom but disre- 
gard interactions. The system is described by two inde- 
pendent Hilbert spaces spanned by the states \0 a ) and 
|1 CT ) with a =t, i- The probability vector for each spin a 
can be written in terms of the eigenvalues and eigenvec- 
tors of the matrix A^ 1 ) (for the two-dimensional Hilbert 
space) as 



-fat 



1 
-1 



(5) 



where the right eigenvector corresponding to the eigen- 
value zero of A^ 1 ) defines the occupation probabili- 
ties for the equilibrium state, P° q = (pg q CT , p cq a ) T = 
(1 - /(e),/(e)) T , with the Fermi function /(e) = [1 + 
exp(/3e)] -1 and the inverse temperature /3 = I/Zcb? 1 . Fur- 
thermore, n a = (0, 1) is the vector representation of the 
number operator for dot electrons with spin er, whose 



initial/equilibrium expectation value is obtained by mul- 
tiplying it from the left into the initial/equilibrium prob- 
ability vector, (rv) in / cc J = n a ■ P^ /eci . The rate 7(T = T 
is obtained as the negative of the non-zero eigenvalue 
of A^ 1 ), with the corresponding left eigenvector being 

l*= (0,l)-(fv) oq (l,l). 

The time evolution of the occupation of each spin state 
is governed by a single decay rate T, 



1 



-rt\ 



,-rt 



(6) 



This equation can be obtained making use of the fact 
that the time evolution of the expectation value of any 
operator, which describes an observable of the QD, is 
given by projecting its vector representation from the 
left onto Eq. ©. 

The time evolution of the total charge of the dot, 
(n)(t) = (fif)(t) + (ni)(i), is also determined by a sin- 
gle relaxation rate j a = T. This means that both charge 
and spin, which are quantities related with single-particle 
processes, do not evolve independently from each other 
and the corresponding decay is given by the same rate. 
A similar non-interacting problem has been studied non- 
pertubatively in Refs. and[H. 

As a next step we consider the squared deviation of 
the charge from its equilibrium value, [n — (n) cq ] 2 . Its 
time evolution is obtained from Eq. (|5|) as 

(ih-(hrf)(t)-m cq } 2 (7) 

= £ [l + (n a r](h a )(t) + 2(n t n i )(t) 

<7=f,4. 

The last, two-particle term of this expression exhibits a 
decay rate given by exp(— 2Tt). This is in contrast to the 
spinless case, where such a term does not appear since 
double occupation is not possible. 

Such an additional exponential decay with the rate 2r 
appears directly in the time evolution of the probabil- 
ity vector, when considering the full two-particle Hilbert 
space spanned by the basis {|0), | f), | 4-), \d)}- In this ba- 
sis, Eq. ((3]) for the non-interacting regime can be written 



/ - [i - m] \ 

|[l-2/(e)] 



-IV 



((hr - (nD 



m) in - (m) oq ) 



(8) 

where as before, P eq defines the state at equilibrium. 
The decaying part of the probability vector can be di- 
vided into three contributions which appear depending 
on how the initial state at to differs from the equilibrium 
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state. Deviations of charge and spin from their equilib- 
rium value relax with the same rate T. The correspond- 
ing expectation values are calculated by multiplying the 
probability vector Eq. (|5J) from the left with the vec- 
tor representation of the operators h = (0,1,1,2) and 
s = (0, 1, — 1, 0) which represent the charge and spin, re- 
spectively, in this two-particle basis. The two left eigen- 
vectors of the matrix A*- 1 ) with the same eigenvalue — T, 
are given by 1„ = h — (n) cq (l, 1, 1, 1) and l s = s. 

The third contribution to the decay of the system into 
the equilibrium comes from the relaxation rate 2r, which 
enters the probability vector in connection with a quan- 
tity rh, defined by the operator in vector notation 

m = (0,/(e),/(e),-l + 2/(e)). (9) 

The left eigenvector of A^ 1 ' with the eigenvalue — 2r 



Again, P cq is the eigenvector of A^ = W^ 1 ' with 
the zero eigenvalue and represents the equilibrium state 
in lowest order in the tunnel coupling (the explicit 
form of the four-dimensional matrix A^, together 
with its entire set of eigenvalues and eigenvectors, is 
given in the Appendix). In the two-particle basis 
{|x>} = {|0),|t) ) U>,|d>}> again n = (0,1,1,2) repre- 
sents the charge operator, and s = (0, 1, —1, 0) represents 
the spin operator. The form of the operator rh is modified 
by the presence of finite Coulomb interaction; the explicit 
form will be discussed later in this sub-section (see Eq. 
(|14l) below). The initial and equilibrium expectation val- 
ues for these operators, entering in the above Eq. (flUl) . 
are obtained as (6) m / cq = 6 ■ p m /°q 5 with 6 = s,n,rh. 
Explicit expressions for (h) cq and (m) eq are shown be- 
low. The negative of the other three eigenvalues of A^ 
directly determine the decay of charge, spin^ and the 
quantity denoted by rh. These decay rates read 

7 „ = T[l + f(e)-f(e + U)] (Ha) 
7s = r[l-/(e) + /(e + [/)] (lib) 
7 m = 217 (11c) 

Notice that due to interaction, the relaxation rates for 
charge and spin (j n and 7 S respectively) differ from each 
other and depend on the level position e, in contrast to 



is given by rh — (m) cq (l, 1,1,1). In contrast to charge 
and spin, the quantity represented by rh does not have a 
straightforward intuitive interpretation, since it depends 
on the quantum dot parameters at t > to and on the 
temperature and chemical potential of the reservoir via 
the Fermi functions. 



2. Finite Coulomb interaction U 

From now on we assume a finite on-site Coulomb re- 
pulsion U on the dot. Analogously to the noninteracting 
case discussed before, from Eq. (|4]) we can write the time- 
dependent probability vector in terms of contributions 
exhibiting different decay times 



(10) 



the non-interacting case. Their dependence on the level 
position is shown in Fig. O In the region for — U < e < 
0, 7„ is enhanced as the charge decays into the twofold 
degenerate state of single-occupation, whereas the spin 
relaxation in first order in T is suppressed, since spin-flip 
processes are not possible. However, the third decay rate, 
7 m , remains fully energy independent as in the case with 
[7 = 0. 

The right eigenvectors occurring in Eq. (1101) each repre- 
sent a change to the steady state density matrix that de- 
cays exponentially with rate 7.; (i = n,s,m). Therefore, 
a system being brought out of equilibrium by a symmet- 
ric deviation between p^ and p^ only, is decaying with a 
rate 7 S . A deviation from equilibrium in which the occu- 
pation of the even sector, po +pd is symmetrically shifted 
from the odd sector, p^ +Pi, is governed solely by the re- 
laxation rate 7m . This right eigenvector is found to play 
an important role also in the low-temperature renormal- 
ization of this models An energy-dependent change in 
the occupation probabilities as prescribed by the second 
vector in Eq. (fit)]) yields a decay of the total charge of 
the system with the rate 7„. The conditions under which 
specific deviations from the equilibrium state should be 
performed in order to obtain a specific decay rate, are 
discussed in the following Section. 



P(t) = P cq 



1- f{e) +f(e + U) 



( -[WW] \ 
l[l-f(e)-f(e + U)] 
f[l-/( e )-/ (e + t/)] 

V f(e + U) J 



7nt n?,\m _ ^,\eq 



n - n> 



/ \ 
1 

2 

V J 



e -7 s */J\in + 



-1 

1 

1 
-1 
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FIG. 2: (Color online) Decay rates 7 m (blue, dashed line), 7„ 
(red, dash-dotted line) and 7„ (green, solid line) in units of V 
as a function of the dot level position e. The temperature is 
k^T = 1.5F and the interaction energy is U = 10F. 



The attribution of these relaxation rates to the charge, 
spin, and rh arises from the independent decay of these 
quantities, due to the explicit form of the left eigenvec- 
tors of A' 1 ). The spin operator coincides with the left 
eigenvector associated to the eigenvalue — j s and since it 
has a vanishing equilibrium value, the time evolution of 
its expectation value is given by 



(0, 1, -1, 0) 



/Po(t) \ 

Pt(t) 

Pi(t) 
\Pd(t)J 



= (s)(t)=e-^(sr. (12) 



Equivalently, the left eigenvector corresponding to the 
eigenvalue — 7„, is h — (n) cq (l, 1, 1, 1). It contains the 
charge operator n and its equilibrium value (n) cq = 
2/(e)/ [1 + /(e) - /(e + U)\. Hence, for the time evolu- 
tion of the charge we find 



(0, 1, 1, 2) 



/ Po(t) \ 

PtW _ (r 

Plit) ^ 

\Pd(*)/ 

(n)(t)-(hr 

e -7„t ^in _ 



\ eq 



The quantity decaying with the rate 7 m alone is related 
to the left eigenvector rh — (m) oq (l, 1, 1, 1), where the 
operator rh is given by 



1 



l-/(e) + /(e + t/) 



( ° 

f(e + U) 
\-l + /(e)+/(e + t/) 



(14) 



Its expectation value follows a time evolution equivalent 
to the one for the charge in Eq. (|13p : (rh)(t) — (fh) cq = 



-7m* 



((' 



i) oq ). Its equilibrium value (rh) 



cq 



FIG. 3: Equilibrium value of the quantity m as a function of 
the dot level position e. The other parameters are: k&T = 
1.5r and U = 10r. 



/(e)/(e + CO/ [1 - /(e) + /(e + U)], plotted in Fig. El is 
- in contrast to spin and charge - not sensitive to the 
regime of single occupation on the quantum dot. Instead, 
it exhibits a feature close to the electron-hole symmetric 
point of the Anderson model, indicating that rh repre- 
sents a quantity which is affected by two-particle effects 
and it decays with a rate that is not modified by the 
Coulomb interaction U. 

Already for the nonintcracting case, we found that 
the rate 2r appears as a consequence of introducing 
two particles in the system, and we considered the de- 
viations from equilibrium charge as a quantity involv- 
ing two-particle processes leading to such a decay rate. 
Also in the case for finite Coulomb interaction, the time- 
dependent mean squared deviations ([n — (n) oq ] 2 )(£) are 
suitable to reveal the relaxation rate 7 m = 2r. Their 
time evolution is obtained by means of Eq. (jTU]) and reads 



([n-(n) cq ] 2 )(i)-[("> 



oqi2 



C-(n){t)-2-(m){t) 



115) 



where in front of the time-dependent charge (n)(t) the 
following coefficient appears: 



C = -2( Pd - Po ) cq + S 



(13) with 



S 



l-f(e)-f(e + U) 
l-f(e)+f(e+Uy 



(16) 



(17) 



The quantity (p<j — po) 

- [1 - /(e) - /(e + U)] I [1 + /(e) - /(e + U)] is the 
difference between the probability of doubly occupied 
and empty dot in equilibrium, which can also be 
related with the occupation of electrons and holes, 
(n e h) = 2(pd — Po)- The behavior of (n e h)° q is shown in 
Fig. |U For e < — U, when the dot is doubly occupied, 
(f l eh) eq = 2; for — U < e < 0, when one electron and 
one hole are present in the system (singly occupied dot), 



6 




FIG. 4: (Color online) Equilibrium electron-hole occupation 
(^eh) cq (red, solid line) and the coefficient S (blue, dashed 
line) as a function of the dot level position e. The other 
parameters are: ksT = 1.5F and U = 10I\ 



,eq 



0; and for e > 0, when the system is completely 



"filled with holes" (empty dot), (n eh ) eq = -2. The 
quantity S is also shown in Fig. [4] (blue dashed line), 
exhibiting a sign change around e = —U/2, the point at 
which the Anderson model is electron-hole symmetric. 
By replacing e4-6-(/,we go from the electron-like to 
the hole-like behavior, finding an inversion in the sign of 
S, S — > —S. The function S therefore indicates whether 
the spectrum of the quantum dot is electron-like or 
hole- like. 

The mean squared deviations of the charge from its 
value at equilibrium is an example for a physical quanti- 
ties showing a decay with j m ; it also includes the charge 
relaxation rate 7„, which is found independently from 
the time evolution of the charge. Equivalently also the 
time-resolved charge variance, ([n - (h)(t)] 2 )(t), or the 
time-resolved spin variance^ ((s) 2 )(i), contain a contri- 
bution decaying with -f m . 



B. Response to an external perturbation 

We now consider in detail which external perturbations 
are necessary in order to induce a decay of the full oc- 
cupation probability vector with one certain relaxation 
rate only, in a controlled way. Furthermore, we address 
the conditions under which a single decay rate can be ex- 
tracted more easily from the occupation of a single state 
by a measurement with a nearby quantum point contact 
(QPC). 

We first address the case of an infinitesimal perturba- 
tion (linear response). A small variation of the gate po- 
tential leads to a decay of the charge governed by the 
charge relaxation rate 7„. Similarly, the infinitesimal 
variation of the Zecman splitting in the dot yields a decay 
with the spin relaxation rate 7 S . In order to obtain a de- 
cay of the state with the rate j m only, it is not sufficient to 



modulate the gate voltage, also the two-particle term in 
the Hamiltonian, Un^n±, needs to be varied. The on-site 
repulsion U could be changed, for example, by tuning the 
carrier density in a nearby two-dimensional electron gas, 
thereby controlling the screening of the electron-electron 
interaction in the dot. From Eq. (|10[) we know that a dy- 
namics given only by 7 m is obtained if the occupation of 
the even states are changed in the same direction, oppo- 
site to that of the single occupied states; this condition is 
fulfilled if infinitesimal variations of the gate, e — > e + e?e, 
and of the interaction, U — > U + dU, obey the relation: 



dU 



1 + 2 cxp(/3e) + cxp(/3[2e + U] 
I + exp(/3e) 



■de. 



(18) 



This expression is represented in terms of field lines in 
Fig. \5\ An infinitesimal change tangential to the field 
line passing through the point corresponding to the initial 
values of e and U leads to a pure decay with j m . 

For parameter variations that are not infinitesimal (be- 
yond linear response), a change only of the gate voltage 
results in a decay of the state with both rates 7„ and 
7 m . From Eq. (|10p we find that a finite variation of the 
energy level and the interaction from an initial condition 
(eo, Uq) to (e = e + Ae, U = Uq + AU) resulting in a 
relaxation containing solely 7„, satisfies the equation 



(3AU = -13 Ae + ln[2 - e _ " Ae ]. 



(19) 



A relaxation given only by the rate 7„ 
relation: 



is found when the 



U = U + - In 



1 _|_ e /3e _|_ e P(e+e +U ) _ e /3(2t + U ) 



(20) 

is fulfilled. For different values of eo and Uo, Eq. (f2"0)) 
produces again the field lines shown in Fig. [SJ Therefore, 
finite variations of the parameters between two points 
lying on the same field line yield a dynamics for the entire 
occupation probabilities vector P governed only by 7 TO . 
Obviously, a generic variation in both e and U which 
does not fulfill the conditions specified by Eqs. (fl"9]) or 
(|2"0"|) exhibits a dynamics of the probabilities with two 
time scales: 7„ and j m . 

In Fig. [5] it is observed that in the region e > —U/2 
the field lines are approximately horizontal, i.e, only the 
interaction U needs to be varied while keeping the level 
position constant in order to see a dynamics of the prob- 
ability governed by j m only. In fact, in this regime the 
QD is predominantly empty and variations of the inter- 
action strength U do not affect the occupation of the dot. 
This is the reason why this variation yields a dynamics 
in which the rate 7„ does not contribute. On the other 
hand, in the region for e < — U/2 in order to avoid that 
the number of particles on the dot changes, which would 
lead to a relaxation with rate 7„, a variation in U needs 
to be accompanied by an opposite variation in e, that is 
Ae = —AC/. The crossover between the two regimes ap- 
pears around the symmetry point of the Anderson model, 
e = -U/2. 



pu 



rt 



FIG. 5: Field lines describing variations of e and U that lead 
to a response of the system with only the rate j m . 



Importantly, it is also possible to read out either the 
rate 7„ or the rate -f m by varying the gate voltage only 
(and, thus, not fulfilling Eqs. JlgJ and ([20])), which is 
easier to realize in an experiment. This can be done 
by measuring an observable that is sensitive to only one 
occupation probability, for instance the probability of the 
quantum dot being empty. Such a time-resolved read- 
out of the probability can be achieved by considering a 
QPC located nearby the system and tuned such that it 
conducts only if the QD is empty4i~— In the simplest 
model of the QPC, which assumes a very fast response, 
the operator corresponding to the current in the QPC is 
given by 



l QPC 



to (1,0, 0,0), 



(21) 



where io is a constant current, given by the character- 
istics of the QPC potential. The expectation value of 
the QPC current is simply (iQpc)(t) = io(po)(t)- In this 
way, the QPC effectively measures the dynamics of the 
occupation probability po . According to Eq. (|10p . a mod- 
ulation of the gate in which the initial value (m) m equals 
the equilibrium value (m) eq leads to a pure decay with j n . 
Instead, for a decay given by j m either the factor (n) ln — 
(n) cq or the factor [1 - /(e)] / [1 - /(e) + /(e + U)\ in 
Eq. (fl0|) has to vanish. 

Results for the QPC current for different variations of 
the level position e while U is kept constant, are shown 
in the logarithmic plot in Fig. [51 For clarity, we also 
subtracted the corresponding current in the long time 
limit, (i) Gq . In particular, for a fixed value of U equal 
to bksT, we find that if the level position is changed 
from eo = lOfceT to e = — lOfc^T, the time evolution 
of po is governed entirely by the rate j m , giving rise to 
the straight, blue-dashed line in Fig. [51 Its slope is given 
by 7m, making it possible to extract this relaxation rate 
from measurements of the current in the QPC. However 
we can obtain a dynamics of po given mainly by the rate 
7n by performing a variation in e from eo = lOfc^T to 
e = 2kpT which results in the red dot-dashed straight 
line in Fig. [5] again, the slope yields the corrcspond- 



FIG. 6: (Color online) Logplot of the current in the QPC 
as a function of the time after a finite variation of e. The 
on-site Coulomb repulsion U is constant and takes the value 
U — 5fcsT. Dashed blue line: e changes from eo = lOfcsT to 
e = — lOfcsT, its slope yields the relaxation rate 7 m . Red dot- 
dashed line: in this case eo = lOfcsT to e = 2fcsT, and the 
slope leads to y n . The black line is obtained if e changes from 
eo = — 13fc_eT to e = 2fcsT, in which both rates y m and f n 
are present. In all cases we have subtracted the corresponding 
value for the current in the long-time limit. 



ing relaxation rate which takes the value 7„ = 1.12 T. 
Finally, we show an example in which variations from 
6q = —ISksT to e = 2kpT (solid black line) produce 
a dynamics of pa which includes two exponential decays 
with rates 7 m and 7„. As a result, the curve exhibits a 
change in the slope, showing that a single rate will not 
be obtained by arbitrary variations of the parameters. 



C. Second-order corrections in the tunnel coupling 

In the previous sections we investigated the relaxation 
rates in first order in the tunnel coupling strength T. 
However, corrections due to higher order tunneling pro- 
cesses appear when the tunnel coupling gets stronger. 
Besides quantitative corrections, this reveals an interest- 
ing new aspect. In second order in the tunnel coupling, 
the matrix included in the exponential decay takes 
the form A< 2 ) = W^+aW^.W* 1 ). The second-order 
corrections to the relaxation rates for charge and spin are 
given by£2- 

7< 2) = <r(e,r,[/)-7„ + <7 r (e,r,t/)7 n 

f(e + U)Wo d + [l-f(e)}W do 

l-/(e) + /(e + t/) 1 > 

7« = a(e,r,[/)^7 s + a r (e,r,[/)7 s + 2W sf . (23) 

These corrections contain renormalization terms as well 
as real cotunneling contributions. On one hand, 
the renormalization terms contain an effect due to 
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the level renormalization e — > e + a(e,T,U), with 

a{e,T,U) = r[<fte + U)-<f>(e)], 0(e) = ^Re* + 

and Sfr(x) is the digamma function. On the other 
hand, the renormalization of the tunnel coupling ap- 
pears, r T[l + crr(e,r, [7)], with a r (e,T,U) = 
-S [Tcj)'{e) + T(j)'(e + U) - ^a(e, T, U)] and where S was 
defined in Eq. (|17|). Real cotunneling contributions are 
manifest in terms of spin flips, W s f, and coherent transi- 
tions changing the particle number on the dot by 2, Wod 
and Wdo- These cotunneling terms read 



W d0 



Ycj>"{e)+Yct>"{e+U)--a l {e,U) 



(24) 



2r 



aJ(2e+U) _ I 



T ( j>'(e)+T<f ) '(e + U)--a(e > U) 



(25) 



erndWoa =exp[/3(2e + lO]Wd . 

The way in which the cotunneling contributions en- 
ter in the respective charge and spin relaxation rates is 
related to the deviation of the state of the QD from equi- 
librium, given by Eq. (|10p in first order in T. As an ex- 
ample we discuss the correction to the charge decay rate, 
second line of Eq. (|22|). There the factor 2 appears due 
to the change in the charge by ±2 in a process bringing 
the dot from zero to double occupation and vice versa4^ 
The fraction with which the transition from zero to dou- 
ble occupation, Wdo, enters the correction to the charge 

(2) 

relaxation rate, 7„ , is given by the deviation from equi- 
librium of P(i) in the direction of po, of the contribution 
which actually decays with j n only. This is the first com- 
ponent of the second vector in Eq. (JTUJ). Equivalcntly, the 
transition from double to zero occupation, Wod, enters 
with the fraction given by the fourth component of the 
same vector, namely by the deviation from equilibrium 
of P(<) in the direction of pd- 

Strikingly, in contrast to the charge and spin relaxation 
rates, 7 m does not get renormalized at all by second order 
tunneling processes: 



£ } = 7^r(e,r,[/) + 

= o , 



l-f(e)-f(e + U) 
l-f(e) + f(e + U) 



(hi 



W d0 ) 
(26) 



The reason for this is that the contribution due to T 
renormalization and those due to coherent processes be- 
tween empty and doubly occupied dot, cancel each other. 
The lack of second order corrections, confirms that this 
relaxation rate is related to a quantity which is not sen- 
sitive to the Coulomb interaction. The fact that correc- 
tions are missing, is also found using a rcnormalization- 
group approach.—. 

Another important aspect of this missing second-order 
correction is that it is due to an exact cancelation of 
the contribution due to virtual second order processes, 
namely the r-renormalization, with real cotunneling con- 
tributions. This is in contrast to, e.g. the conductance, 



where only the real cotunneling processes contribute far 
from resonances, while renormalization terms are limited 
to the resonant regions. 



D. Hybrid systems 

Until now, we considered the quantum dot to be cou- 
pled to a normal conducting lead. However, the vicin- 
ity of a superconducting or a ferromagnetic reservoir in- 
duces correlations between electrons and holes or between 
charge and spin, respectively. In the following we study, 
in first order in the tunnel coupling strength T, the influ- 
ence of induced correlations on the relaxation rates of the 
dot. The charge response of a noninteracting mesoscopic 
scattering region coupled to both normal and supercon- 
ducting leads has been studied in Refs. lisHitl 



1. Proximity to a superconducting lead 

In the previous sections we have seen that the rate 
7 m , which together with the time decay of charge and 
spin determines the relaxation of the QD to the equilib- 
rium state, is independent of the level position and the 
Coulomb interaction and that it enters in the time evo- 
lution of quantities sensitive to two-particle effects. It is 
therefore expected that the rate 7 m will directly influence 
the relaxation of the charge towards the equilibrium in a 
setup that naturally mixes the empty and doubly occu- 
pied states of the dot. This situation is obtained if the 
QD is not only coupled to a normal lead (with tunnel 
coupling strength T) but also to an additional supercon- 
ducting contact (with tunnel coupling strength Tg), as 
shown in Fig. [T] (b). We consider only the case when the 
superconductor is kept at the same chemical potential as 
the normal lead and we set both chemical potentials to 
zero. The only purpose of the extra lead is here to induce 
superconducting correlations on the dot via the proximity 
effect. To the original Hamiltonian, Hr> + H rcs + Ht, we 
now add the Hamiltonian for the superconducting con- 
tact and its tunnel coupling to the QD, 



H s = ^e S fe4 feff c Si . .-^(Ac s _ fc4 _c Si . t + H.c.) 

ka k 

E^4a^+H.c). (27) 



A;, it 



where c^l a is the the annihilation (creation) operator of 
electrons in the lead. In the limit of a large supercon- 
ducting gap A the effect of the additional contact can 
be cast in an effective Hamiltonian of the dot which in- 
cludes a coupling between electrons and holes in the QD, 
H [ ° S) =H D - T s /2(d\d\ + H.c). The eigenstates of the 
proximized dot are the states of single occupation |cr) and 
other two states which are superpositions of the empty 
and double occupied states of the dot (due to Andrccv 
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reflection) 



(28) 



with energies given by i?± = 5/2 ± e a, where the level 
detuning between |0) and \d) is 5 = 2e + U and 2ea = 
yPTTf is the energy splitting between the |+) and 
|-) states In the new basis {|+), | f), 1 1), \~)}, the 
vector representing the charge operator is expressed as 

h = ^— .^2 + j^, 1, 1, ^2 — -^j and we expect that the 

effect of the mixing of electrons and holes will be visible 
in its time evolution. In first order in the tunnel-coupling 
strength to the normal reservoir Y and assuming r <C Ts, 
we find the relaxation rates 



7s,i = T[l 

75,. = r[i 
is,2 = 2r. 



f(e- 



-E-) 



-f(E+ 
f(E+ 



0] 



(29) 
(30) 
(31) 



Remarkably the eigenvalue — 2r = —75,2 remains un- 
affected, i.e. 7„i = 7s,2 is not modified by the presence 
of the additional superconducting lead. 

The spin on the dot, which is determined by the oc- 
cupation probabilities of singly occupied states, still de- 
cays with a single relaxation rate given by "fs,s, i-e. 
l s = (0,1,-1,0) is an eigenvector of the kernel A^ 1 ) 
(in the proximized basis). In contrast, the decay of the 
charge to its equilibrium value is given by 

(h)sc(t) = \ [{n) in - (n) cq ] (e"^ 2t + e -wt) + ( - )oq 



-a sc - [(x) in - (x) 



+2 [(y)' m - <y> oq ] ( e_7s ' 2t - e_7s ' 14 ) ( 32 ) 



with 



asc 



(2 



f(E- 



'f(E + -e)+f(E_-e) 



and where we defined the difference in the occupation 
of the |±) states, x = p+ — p- and the quantity y = 

{k- — T)p+ + (fc+ — l)p~, with k± 



T 



2 ±-2-. The 



charge evolves with two different time scales, 75^ and 
7s,2 = 7m, instead of only one as in the normal case. This 
is a direct consequence of the mixing of the states |0) and 
\d) induced by the superconducting contact. This effect 
opens the possibility to extract this rate by measuring 
the time evolution of the charge in the proximized dot. 



not been modified. Since we associate this rate with pro- 
cesses involving two particles each with spin tr, it is ex- 
pected that if the spin symmetry is broken by introducing 
a ferromagnetic contact, the rate 7,„ will now be the sum 
of the tunneling rates for spin up and spin down electrons. 
In order to verify this, we consider the Hamiltonian used 
for the normal case and assume a spin-dependent density 
of states in the only reservoir attached to the quantum 
dot, sec Fig. [T](c). This leads to spin-dependent tunnel 
couplings, I 1 -)- and T^, which are included in the corre- 
sponding transition matrix . Diagonalization of A^ 
yields the three relaxation rates: 

7F.i = r + y (AT) 2 + 41>r; [/(e) - fje + U)p3) 



if.2 = r - y (at) 2 + 4r t r ; [/(e) - f(e + U)p4) 



2T 



(35) 



with r = § (T i + T t ) and Ar = I> - T t . 

As in the normal case, there is an eigenvalue which does 
not depend on the level position nor on the interaction 
but on the sum of the different tunneling rates: — 2r — > 
— (r^ + I 1 ^) . The appearance of such a combination of 
the spin-dependent tunneling strengths in the relaxation 
rate, confirms the statement that two-particle processes 
involving electrons with both spin polarizations are at 
the basis of the decay rate 7 m . 

Due to the ferromagnetic lead, the dynamics of spin 
and charge are now mixed. The corresponding time evo- 
lution in first order in the tunnel coupling takes the form: 



+a s (s) ln (e 
+b s [(f 



3\in/ -7jr,it 



-7-F,2*"j 



(36) 



-7f,2*1 



1 



(n) F (t) = - [{h) in - (n) cq ] (e"W + e^ 2 ') + (n) cq 
+a c [(n) in - (7i) cq ] (e-w* - e^ 2 *) 
+Ms) in (e- W -e- 7F < 2t ) (37) 
where we introduced the abbreviations: 



b s = 



b c = 





r[/(e)- 




CO] 






/Ar 2 + 4r t r^ 


im 


-f(e + 


CO] 2 




Ar[l + /(e 


>-/(■ 


z + U)] 




2 V 


/ Ar 2 + 4r T r i 


[/(e) 


-/(e + 


U)} 2 


a s 


Ar[l - /(e 


> + /(' 


e+CO] 




2 Y 


/Ar 2 + 4r t r 4 


[/(e) 


-f(e + 


U)} 2 



2. Ferromagnetic lead 

Even though the presence of a superconducting lead 
couples electrons and holes, the relaxation rate "f m has 



The last term in Eq. ([31))) shows that at finite time t the 
initial charge influences the time evolution of the spin; 
similarly, the initial spin enters explicitly in the dynamics 
of the charge, Eq. (|3"T|) . These terms persist in the non- 
interacting limit, revealing that the coupled evolution of 
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charge and spin including two relaxation rates (which for 
the non- interacting case take the form r fF\ = and 
7i?2 = TjJ is a direct consequence of the presence of the 
ferromagnetic contact. In contrast, the factor a s = a c 
vanishes for U = implying that it stems from the com- 
bined effect of the Coulomb interaction and the break- 
ing of the spin symmetry. As expected the independent 
evolution of charge and spin is recovered in the limit 
= r^. The mixing of the dynamics of both, charge 
and spin, induced here by a ferromagnetic lead was found 
in Ref. [29| for the case of lifted spin-degeneracy in the 
dot due to a finite Zeeman splitting. Note that for the 
hybrid as well as for the normal system, the sum of the 
energy-dependent relaxation rates equals 21" , as long as 
the tunnel coupling T is treated in first order, only. 

IV. CONCLUSION 

We have studied the different time scales present in 
the evolution of the reduced density matrix of a single- 
level QD with Coulomb interaction and tunnel coupled 
to a single reservoir, after being brought out of equilib- 
rium. Besides the relaxation rates for charge and spin, 
we find an additional rate "f m = 2T, which is independent 
of the energy level of the dot as well as of the interaction 
strength. This relaxation is related to the presence of 
two particles in the dot and is found to be not sensitive 
to the Coulomb interaction. The time evolution of the 
square deviations of the charge from its equilibrium value 
is proposed as a physical quantity related with processes 
involving two-particles leading to the rate 2T. In order 
to further elucidate the properties of this decay, we ana- 



lyzed the response of the system to specific variations of 
both, the interaction strength U and the level position 
e, finding that 7 m can be extracted from time-resolved 
measurements of the current passing through a nearby 
quantum point contact. Additionally, we analyzed two 
other setups: a dot proximized by a superconductor and 
coupled to a normal reservoir, and a dot coupled to a fer- 
romagnetic lead. In the hybrid normal-superconducting 
systems, we found that the time-resolved read-out of the 
charge represents another possibility to get access to the 
rate j m . 
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Appendix: Normal case. Eigenvalues and 

eigenvectors in first order in the tunnel coupling 

The transition matrix for the normal case in the eigen- 
basis of the isolated QD {|0), | t), I I), \d)}, in first order 
in the tunneling strength T, is calculated by means of 
Fermi's Golden rule and is given by: 



A d) = w (1) (z = o) = r 



-2/(6) l-/(e) 1-/(6) \ 

/(e) -[l-/(e) + /(e + C/)] l-f(e + U)) 

/(e) _[i_/( e ) + /( e + [/)] l-f(e + U)) 

f(e + U) f(e + U) _ 2 [l-/(e + [/)]/ 



(38) 



with the Fermi function f(x) = 1/[1 + cxp(/3x)], where (3 is the inverse temperature /3 = (fcgT) . 

As A*- 1 ) is non-Hcrmitian it has different right and left eigenvectors, and U. For a system with a well-defined 
steady state (as the one we are considering here) there must be at least a zero eigenvalue, Ao = 0J£ The other 
eigenvalues are found to be the negative of 



7n = T [1 + /(e) - /(e + U)] 

ls = T[l-f(e) + f(e + U)} (39) 

7m = 2r. 



The right eigenvector corresponding with the zero eigenvalue, i"o, determines the stationary density matrix (which we 
also label as P cq ), whereas each one of the rest of the right eigenvectors represents a deviation out of the equilibrium 
density matrix which decays exponentially with a rate given by the negative of the corresponding eigenvalue: 
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ro = 



[l-/(e)][l-/(e + C/)] 
f(e)[l-f(e + U)] 



1 + f(e)-f(e + U)\ /(e)[l-/(e + f/)] 

/(e)/(e + EO 



/ 



l-f(e) + f(e + U) 



-[i-mi \ 



l[l-f(e)-f(e + U)} 



i 



2 [l-f(e)-f(e+U)] 
f(e+U) 




(40) 



The vector of equilibrium occupations ro can equivalently be written in terms of the Gibbs factors as 

1 



ro = P cq = 



/ 1 \ 

P -Pe 



1 + 2e-P* + e-Pi^+u) 



\ 



-/3(2e+C/) 



(41) 



The left eigenvectors determine the quantities decaying into equilibrium with a single time scale only, and are found 
to be 



lo = (1,1,1,1) 

l s = (0,1,-1,0) 

1„ = (n)«» (1,1, 1,1) -(0,1, 1,2) 



1, 



(mr (1,1, 1,1) 



l-f(e) + f(e + U) 



(0, /(e + U), f(e + U), -1 + f(e) + f(e + U)) 



(42a) 
(42b) 
(42c) 

(42d) 



where cc " = 2/(e)/[l + /(e) - f(e + U)} and (to) cc i = 
[f(e)f(e + U)]/[l-f(e) + f(e + U)]. 

Normalization constants are chosen such that the 
eigenvectors fulfill the relation lj • = (S^r^ with i, j <E 
{0, s, n, to}. 

These left eigenvectors contain the operators for spin, 
charge and to in vector representation, which can be 
understood in the following manner. While in general 
the expectation value of an operator O is found from 

(0)(t) = Tr |Op(i) j, with the full density matrix p, this 

can be considerably simplified in the situation considered 



here, where only diagonal elements of the reduced den- 
sity matrix of the quantum dot, collected in the vector 
P, play a role. The expectation value of a quantum dot 
operator is then obtained by multiplying its vector rep- 
resentation from the left hand side onto the vector P. To 
show an example the expectation value of the spin on the 
dot is obtained by multiplying P from left by the vector 
(0, 1,-1, 0), yielding (s) = 1 ■ p t + (-1) ■ p i . Similarly, 
all other operators for quantum dot observables can be 
expressed in such a vector representation. 
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